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Abstract 



We study spherically symmetric solutions in f{R) theories and its compati- 
bility with local tests of gravity. We start by clarifying the range of validity of 
the weak field expansion and show that for many models proposed to address the 
Dark Energy problem this expansion breaks down in realistic situations. This 
invalidates the conclusions of several papers that make inappropriate use of this 
expansion. For the stable models that modify gravity only at small curvatures we 
find that when the asymptotic background curvature is large we approximately 
recover the solutions of Einstein gravity through the so-called Chameleon mecha- 
nism, as a result of the non-linear dynamics of the extra scalar degree of freedom 
contained in the metric. In these models one would observe a transition from 
Einstein to scalar-tensor gravity as the Universe expands and the background 
curvature diminishes. Assuming an adiabatic evolution we estimate the redshift 
at which this transition would take place for a source with given mass and radius. 
We also show that models of dynamical Dark Energy claimed to be compatible 
with tests of gravity because the mass of the scalar is large in vacuum {e.g. those 
that also include corrections in the action), are not viable. 
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1 Introduction 



Theories of gravity whose action is some function of the Ricci scalar have received much 
attention recently as possible models of dynamical Dark Energy. It is well known that 



to find R{4>) for a range of R, a solution where R lies in this range can also be found 
from an equivalent action that consists of Einstein gravity minimally coupled to a 
scalar field, namely 0, with certain potential [1|2J. So these models, when applied 
to the Dark Energy problem, are equivalent to quintessence models where the scalar 
field is co,*n„a.,v coupled to .natteg It has been shown that for many choices 
of the function / one obtains acceptable cosmological solutions that could be used 
to describe the acceleration of the universe: in fact one can reproduce an arbitrary 
expansion history by choosing / suitably [4.5]. What has been however subject of 
controversy is whether one can, at the same time that one describes the Dark Energy 
in this way, satisfy the bounds on the existence of extra scalar fields with gravitational 
couplings coming from Solar System or laboratory experiments (see e.g. [2|6|7f8] ). 

Let us explain here the reasons why this controversy has arisen. In order to check 
the experimental implications of a gravitational theory for weak fields one usually 
starts by considering a background solution and making an expansion of the equations 
of motion (EOM) in powers of the fiuctuations over this background solution. If the 
second (and higher) order terms in this expansion are negligible with respect to the 
first order ones in a considered solution, perturbation theory is applicable. When this 
is the case, the EOM can be approximated by a set of linear second order differential 
equations, and we can usually find analytic solutions for static and symmetric sources 
and check the compatibility with experiments. The situation is however completely 
different when this expansion is not applicable. In this case we have to deal with non- 
linear differential equations that are usually difficult to solve, and addressing the issue 
of compatibility with experiments becomes challenging. We will see that this is indeed 
the situation for many of the considered f{R) gravities: perturbation theory breaks 

^In this paper we restrict ourselves to the conventional metric formulation of f{R) gravity and do 
not consider the Palatini formulation [3]. 



when we take a gravitational Lagrangian as Cg = Mpf{R)/2 (where Mp is the Planck 
mass and R the scalar curvature), as long as one can invert the relation 
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down in some cases and is not applicable because non-linearities in the equations are 
not negligible. 

In |2|6|f7] this linearized expansion was considered for some f{R) theories. As we 
said this theory generically consists of Einstein gravity plus an extra scalar field confor- 
mally coupled to matter, so at the linear level the only terms that appear in the EOM 
for the scalar are the kinetic and mass terms with a source given by the trace of the 
energy-momentum tensor (EMT). It was therefore argued in [2|l6|f7] that simple models 
in which the mass of the scalar is very small in vacuum, such as / = i? ± fi"^"^"^ / R"' , 
are ruled out because it would mediate an extra long range force contradicting exper- 
iments^. But this analysis was challenged in [8] where solutions with R ^ fi^ where 
considered. Intuitively, we can expect by looking at the action above that when R^ fi^ 
the effect of the modification should be negligible. The extra terms in the equations 
will be suppressed by powers of the ratio fi'^/R and the solution for any given source 
should be very close to that of General Relativity (GR) up to these small corrections. 
But while this is true, we will see in this paper that in fact there is no contradiction 
with the previous analysis. This is so because, in the models with inverse powers of 
R, whenever the scalar curvature deviates from its vacuum value significantly in some 
region {AR/Ro > 1 with Rq ~ fi"^), the linearized expansion breaks down, and the 
analysis of [2|6|7] can not be applied to those solutions. It is then conceivable that 
in these models an effect similar to the Chameleon mechanism of [9] could take place, 
where the effects of a scalar field that is very light in vacuum are hidden because 
non-linearities are important when considering sources. As we will see, the asymptotic 
boundary conditions for the scalar curvature are the crucial element that will allow 
us to decide whether there is a Chameleon effect or not for the extra scalar degree of 
freedom. In particular, for the models with inverse powers of R, if the background 
curvature goes to a small value asymptotically {Rq ~ /i^), the Chameleon effect does 
not take place: astrophysical bodies do not have a "thin shell" for static solutions. 
In this case the perturbative solutions of pp|7] are valid and the scalar curvature is 
"locked" to the background value everywhere, ~ -Rq ~ /^^ (a fact also noticed in 
[7|). However, when the asymptotic value of the curvature is Rq ^ fi^, the Chameleon 
effect does take place and we recover the solutions of Einstein gravity up to small cor- 



^In these cases the mass of the scalar in vacuum would be ~ /i, and the application of these theories 
to the Dark energy problem demands /x ~ iJo, where Hq is the value of the Hubble constant today. 
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rectionci- We get then the following picture of gravitational dynamics in these theories 
when applied to the Universe: one would observe a transition from Einstein gravity at 
early times (when ^ /i^) to scalar-tensor gravity at late times (when ~ /i^). 
If we assume that this transition is adiabatic, i.e. the solution is always taken to be 
the equilibrium one, these types of modification are ruled out, since we would be in 
the scalar-tensor regime at present. However, in this paper we do no not rule out the 
possibility that there are some models for which this transition is non-adiabatic and 
slow enough so we would still remain in the "GR regime" locally in the Solar System. If 
this was the case those models would not be in conflict with local tests of gravity. The 
Chameleon effect will also allow us to resolve the apparent discontinuity in the General 
Relativistic /i — ^> limit of these theories where we seem to end up with a massless 
scalar field coupled to gravity instead of GR. In the real Universe {i.e with a non-zero 
background energy density) we would recover the usual gravitational dynamics in this 
limit only at a non-perturbative level. 

On the other hand, it has been shown that for some specific forms of the action 
the scalar field can be made arbitrarily massive in vacuum [T11T2] . If this mass is large 
enough it is claimed that the scalar would have passed undetected so the models do 
not conflict with local tests of gravity. But as we will see this conclusion is wrong and 
this possibility can not be realized. The reason for this is that we can not trust a weak 
field expansion in this case: as we raise the mass of the scalar field, we lower the energy 
scale where non-linearities become important. We will find that for most functions / 



inearities become 
If this mass is 



that attempt to describe Dark Energy the energy scale at which non- 
relevant is As ^ Mp {^f^ ; where nxg is the mass of the scalar fielc^ 
large compared with Hq, as it should be in order to avoid conflict with experiment, this 
scale is very small which means that in all realistic situations where we want to apply 
the linearization over vacuum it is not possible to do so because non-linearities in the 
equations are not negligible. So in this case the fact that the mass term appearing in 
the linearized EOM for the scalar in vacuum is large does not imply compatibility with 
experiments. And we will argue that, for observationally relevant distances, the true 



^We are assuming here that the mass squared of the scalar is positive in this regime, so the solutions 
with i? ^ are stable. For instance, when f = R± this would be the case only for the 

positive sign, see [10]. We will come back to this issue in section four. 

"^And in all models that satisfy some "minimal requirements" to be defined precisely later, it is at 

least of order ~ Mj, 
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solution on vacuum is again of the scalar-tensor type with an effectively massless scalar, 
as in the generic models. Furthermore, we will show that models that include positive 
powers of R to raise the scalar mass in vacuum will also conflict with experimental 
results on large curvature backgrounds, where R ^ /i^. In those backgrounds the mass 
of the scalar diminishes to a small value and the strong coupling scale raises so 
one can use the linearized equations and it is easy to conclude that one should have 
observed the effects of the scalar on laboratory searches of fifth forces, for instance. 

The paper is organized as follows. In the next section we will review the solutions 
of the linearized EOM for a spherically symmetric mass in f{R) gravity and we will 
discuss the range of validity of the weak field expansion for a generic action f{R). In 
the third section we will focus on theories designed to explain the cosmic acceleration 
and we will consider the expansion over vacuum where i?o ~ /^^ ~ Hq. We will 
distinguish two cases: the "generic" one (in which the scalar mass is ~ /j.) and the 
"fine-tuned" one (in which the mass is ^ /x) . In the first case we will comment on the 
curious property that the scalar curvature remains small in the perturbative solutions 
even in those regions where there is a large energy-momentum density. In this case 
the strong-coupling scale of the linearization is the Planck mass, so once the curvature 
gets locked into this phase it remains there: one would require Planck scale energies or 
very strong gravitational fields in order to excite the scalar curvature out of its vacuum 
value. In the second case we will explain why when we raise the mass of the scalar in 
a certain background in a model of dynamical Dark Energy, the linearized expansion 
over such background is actually of no use. And we will argue that the true solution 
exhibits the same behavior as for the generic case. Furthermore we will explain why 
the addition of positive powers of R to the action to yield this effect in vacuum would 
imply a conflict with experiment whenever R^ /i^. In the fourth section we will briefly 
discuss the Chameleon mechanism and the recovery of Einstein gravity in models that 
modify gravity only at small curvatures ~ /i^, as long as the asymptotic background 
curvature is Rq ^ //^ and the scalar mass squared is positive. We will also discuss the 
transition from GR to scalar-tensor gravity (or "locked curvature" phase) that would 
take place in these theories when applied to the Universe as the background curvature 
diminishes to a value ~ //^ and the Chameleon effect gradually disappears. Finally we 
end with the conclusions in section five. 
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2 The short distance weak field expansion, general 
case 



It is well known that the EOM derived from the action 

S = / d*xy^f{R) + S-^ (with Gn = 1/(8^)) , (2) 

where Sm is the matter action minimally coupled to the metric, are equivalent to those 
of Einstein gravity minimally coupled to a scalar field that is conformally coupled to 
matter. This equivalence is most easily demonstrated directly at the level of the action 
[Tp] . but it is also instructive to take a slight detour and show how the equivalence 
arises at the level of the EOM. The EOM derived from ([2D read: 

E^^g^u) = f\R)R,. - f{R)^-f + 9^.unf\R) - V,VJ'{R) = StcGnT,, • (3) 

A striking difference with Einstein gravity is seen from the trace: 

Ejiig^,) = -2f{R) + f\R)R + 3n/'(i?) = SttG^T . (4) 

In contrast with GR, where the Ricci scalar is completely determined by the trace of 
the energy- momentum tensor {R = —SttGnT) we now see that R, or rather f'{R), 
is an independent propagating degree of freedom. It will therefore be convenient to 
introduce a new variable ip, that we will identify with f'{R), and cast the original set 
of fourth order differential equations in g^^, into a set of second order equations in g^^, 
and ip: 

E^,{g^,,tP) = ^i^^,-/(7^(^))^ + ^^,□V;-V^V.V^ = 87rGivT^,, (5) 
W) = R, (6) 

where TZ is defined as the inverse of /': 

f{n{^|J))^^|J. (7) 

It is clear that a solution of eqs. flHfB]) also solves the original EOM as long as 7l{ip) 
is well defined on the solution. Notice that the invertibility of /' is equivalent to the 
more familiar condition that f"{R) should be different from zero. In general f"{R) can 
become zero for some values of i? so 7^ will be a multivalued function, and one has to 
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chose a certain branch. The solutions of ([3]) are found then by considering eqs. fl5]|6l) 
with all possible branches. 

Finally, it turns out to be convenient to define a new metric variable 

9^lv = f'{R)9iiv = -ipgiiu , (8) 

since in terms of this metric the eqs.( !5]|6i) let themselves reshuffle as those of Einstein 
gravity, minimally coupled to the scalar ip and with a metric g^i, = g^^/ip in the matter 
sector: 

E,, = ij-'E^,{g,,,ij) + ^{n{ij)-R) 



G 



Svr^T^., (9) 



= E;^(^7^„V^)+V^(7^(V^)-i^)=3□^^ + Vp(^^) = 87^G^^T, (10) 
where (5^,^ is the Einstein tensor corresponding to g^u and we have defined 

Vpiij) ^ ^7^(^) - 2/(7^(^)) . (11) 

One can now verify that these equations indeed arise from the variation with respect 
to §1^1, and ip of the action 

Se, ^ I d'x^g {k - A^d.^jd^ij^ - v(^)) + sU9,um , 

= j d'x^g i-^R - -d,<pd,rg^^ - V{<P)\ + S^{g,,e-^'f>) , (12) 

where 

VW = i(KW-ffl)). (13) 

and in the last line we have reparameterized the scalar: 

0= y|Mpln^ = /€-Mn7/>, (14) 



to get a Lagrangian with the usual kinetic term and a potential V(0) = MpV{e'^'^) /2. 

We are interested in the situation where there is some local fiuctuation on a back- 
ground characterized by some curvature Rq and a density Tq. By local we mean that 



the distance scale of the fluctuation is short with respect to the characteristic distance 
scale associated with the curvature of the background, so that we can effectively take 
the flat space limit. However, as we will see shortly, this does not mean at all that we 
can forget about the background curvature i?o- In the end, the reason for this is that 
the local fluctuations have to match asymptotic boundary conditions that obviously 
do depend on the background. 

The weak field expansion now consists of Taylor expanding the eqs.f pllTOl) in powers 
of the fluctuations and solving them order by order. So we write: 

^ = ^o(l + V^), 

9,iu = Mglu + V) = Mvf^iy + V) , (15) 

where the sub/superscript stands for the background value, and as we said, we take 
the flat space limit (7°^ ~ rj^y. At first order, the scalar equation flTU]) now only depends 
on the scalar fluctuation 'ijj and the trace of the fluctuation of the EMT T = rj^^T^y 

{dl-V^)^ + ml^ = -^-i,Gff, (16) 

where 

2_ 1,,,,,,. ^ ^(f'(R<s) 



™i--3V;W«) = 3^^-fl„j, (17) 

and we have defined the effective Newton's constant G]^^ = Gn/iPq. Later we will also 
use the effective Planck mass defined as: M^^^ = tjjoM^. 

Given the matter source T one can now unambiguously determine ip at linear order. 
So let us review here the case of a spherically symmetric matter distribution with con- 
stant density T = —Too = ~P, extending over a radius r©. The static and spherically 
symmetric general solution of ([T6|) outside and inside the distribution reads: 



i^out = Ci—^ + C2-^, (18) 
i^.n = C,'- + C4^ + ^^!f^. (19) 

One then determines the constants Ci by imposing the appropriate boundary condi- 
tions. First of all we have the condition that ip takes the background value at infinity: 



^Since we are only considering local fluctuations, we can drop terms like r° /i^"^ or d^h'^'^d^ipo- 



^/'(r) — > as r ^ oo. This sets C2 to zero. One then also has the condition that 
ip is regular in the origin, this sets C3 = — C4. Finally one can determine the other 
two constants by matching the solution outside the distribution to the solution inside: 

'ipoutird,) = ipiniro) and ipoutif&) = i'Lif®)- In the case that < r^^, one finds in 
this way that: 

^out - '±J^e-^\ (20) 

where M = AttpTq/S is the total mass of the distribution. Notice that when we can 
neglect the scalar mass term in (JT6II . even for more general static spherically symmetric 
density distributions one can use Gauss's Law to find the same solution outside the 
source with M being indeed the total mass of the distribution. When we can not neglect 
the scalar mass the solution outside will depend on the actual density distribution 
inside. For instance, in the constant density case, the matching of the solution outside 
to the one inside gives for ^ Tq^: 

^out - l!L^!e^..(^o-), (22) 
SrUg r 

-^^^ .»„..- 4^ . _ re^„.„-,„^) , (23) 



3m^ V r 

We are assuming here that the background is stable so that ml > 0. (In fact ~ —Rq 
would also be fine since for our purpose, the study of local fiuctuations, this would mean 
rUs ~ effectively.) 

We have now found the first order solution for the extra scalar degree of freedom -0 
contained in the metric. The first order solution for the fiuctuations h^^^ of cjf^i, follows 
from the linearization of (jH]) which gives the same equations as those of Einstein gravity 
(with a rescaled Newton's constant): 

G2 = 87rG''J%u . (24) 
So for a static spherically symmetric mass distribution we find the usual result: 

^00 ^ —^^ , (25) 

^ 5.,^^^, (26) 



. (27) 
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We can now finally write down the linearized result for the actual metric 

9,iu = ^ ^ Vfiiy + - i'Viiu- (28) 

And in the small scalar mass case, -C r^^, we find for the metric outside the 
distribution the familiar result of a Brans-Dicke scalar-tensor theory with vanishing 
Brans-Dicke parameter u: 

ds^ = dx^dx'^gf,^ (29) 

- ?^ - ?^e-».A + (l + ?^ - ,1.' . 

r 3r / \ r 3r / 

However, this first order result will only be a good approximation if the weak field 
expansion makes sense. That is, if the neglected higher order terms in the equa- 
tions are subdominant to the individual first order terms that we have solved for. 
The weak field expansion of the kinetic terms in the equations (MTOj) is the same 
as for GR, so we already know that this expansion will only break down as Tq ap- 
proaches the Schwarzschild radius ~ G^-^M, where h^^,ijj ~ 1. The additional 
ingredient now is the expansion of the potential term Vp{ip){= —ip^V'{4j)) in the 
scalar equation (ITO!) . Comparing individually the first order terms that we are keeping 
ipQif)" {r) , ipQip' {r) / r , ipQ-m^iljlr) with the higher order terms coming from the Taylor ex- 
pansion of Vp{ip) that we are neglecting, we see that the weak field expansion will be 
valid as long as l 

^^^^(l + m2r2)»(V;o^(r))"V(")U=^, for every n > 1 . (30) 

For all the situations that we consider in the next sections, the terms coming from the 
second and third derivative of Vp will pose no real limitation on the applicability of the 
weak field expansion, and the breakdown of the linearization will be due to the higher 
order terms, that typically become important at short distances or high energies. It 
is then customary to describe the breakdown of the weak field expansion in terms of 
a strong coupling scale A,,. One therefore considers fiuctuations characterized by a 
single energy scale E, such that h^^,,!/) ~ E/M^^^ and d'^h^u,d'^ijj ~ E'^/M^^L The 
strong coupling scale is then defined as that energy scale where the weak field expansion 



^Throughout this paper, when writing expressions hke A ^ i?, we will always mean \A\ ^ \B\ 
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breaks dowiu. From GR we know that the expansion of the kinetic terms will be valid 
as long as E <^ M^^f . While from the expansion of the potential term in (ITOl) we now 
find the strong coupling scale 



A, = min< M^-^, 

n>3 



/ i- J- \ 11—1 

1 / M;ff 



(31) 



Notice that this strong coupling scale would be the same as we would find by simply 
Taylor expanding the potential of the canonically normalized field defined in (JT2l) over 
the background value 0o(= «^^^ln'?/'o) as 

v{ct>) ^ rhs^' + X4>' + X4>' + Y1 wy- (32) 

~ 1 /2 

where (j) = cj) — 0o- The strong coupling scale is then = ipQ x min{Mp, |A.„|}, where 

1 /2 

the factor ^/'q converts the scales in the so called Einstein frame, described by the 
metric g^y ~ i'orj^u, to the corresponding scales in the matter frame, described by the 
actual metric g^j,^ ~ r]^^. 

3 The spherically symmetric solutions on low cur- 
vature backgrounds 

In this section we will focus our attention on f{R) models that intend to explain the 
cosmic acceleration through a modification of gravity at low curvatures. In particular 
we will apply the results of the previous section to check under which conditions we can 
use the linearization of the action over a background of small curvature {Rq ~ H^) to 
compute the solution that corresponds to a spherically symmetric mass distribution. As 
we said we will impose some conditions on f{R) so that one can talk of the acceleration 
as a "curvature effect", as opposed to e.g. the ACDM model that could be described 
by / = -R — A. In particular we will assume that f{R) = R + F{R) where there is 
a single mass scale appearing in F that we will denote by fi and is of order ^ Hq. 
Then we will consider separately two cases: the generic and the fine tuned one. In the 
generic case we will assume that F is such that 

rp(n)\ ,2—271 tt2— 2n (nn\ 



^This scale is also relevant when computing for instance quantum corrections, since for momenta 
higher than this the loop expansion will also break down. 
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for all n > 1. For n = 1 we mean by this that -F'lij^/^a ~ /'|^^^2 ~ 1. For instance, this 
is the case of modifications to GR where F = fi'^"'^'^ / or F = ^^Log{R). We will see 
that in this case we can use the linearization in vacuum to study the metric in all the 
conventional weak field situations of GR, and the solutions conflict with observations. 
In the second case, the fine-tuned one, we will consider the particular situation in which 

F"\r=i2hI « /i'' ~ (34) 
but still F{R) satisfies (15^ for some n > 2. This case includes models like [TT] 

^(^) = -^ + «^, (35) 



F{R) = ^?Log{R) + a— and F{R) = + a—, (36) 



or like ^ ^ 

with a an order one parameter, that have been claimed to be compatible with local 
tests of gravity. In this case the linearized expansion over the background with R = 
12Hq breaks down as we approach any source and therefore can not be used to find 
the solution. We will nevertheless see that the true solutions for these models are 
essentially the same as in the generic case, with an extra massless scalar in conflict 
with observations. 

3.1 Generic case 

In the previous section we learned that the mass of the scalar and the range of validity 

of the weak field expansion follows from the derivatives of Vp{ip) on the background. If 

the condition fl33l) holds one can easily check that ipQ will be an order one number and 

that Vp{ipo) and all its derivatives are of order /i^ ~ Hq. We then find that the scalar 

is essentially massless ~ Hq and the linearized solution fl2I?]) corresponding to a 

static spherical mass distribution is the one of a massless Brans-Dicke scalar field with 

vanishing uj parameter, clearly ruled out by Solar System tests {uj > 40.000 according 

to the latest measurements from the Cassini mission [13j). However, as we said, one 

should still check if the linearization is applicable, since the weak field conditions are 

not necessarily the same as those of GR. But in this case we find that they actually 

are: from fl5U]) we find that the linearization is valid as long as the usual condition 

GnM , , 

^— < 1 (37) 
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is satisfied. So we can use the weak field expansion for the same situations as in GR 
and in those cases the solution (1291) does indeed describe to a very high accuracy the 
static solution that corresponds to a spherical mass distribution and matches the low 
curvature background at infinity. In fact, we see that when translated into an energy 
scale, the effects of non-linearities would not be important until we reach the Planck 
scale since ~ Mp , as happens in conventional GR. There are however significant 
differences with respect to GR. Perhaps one of the most striking features of these 
solutions is that the scalar curvature remains of order i? ~ /i^ even inside sources 
where GnT ^ yU^, as also noted in [7j. But this is easily understood by going to the 
definition of the field ip. For instance, in the case when F = /x^""''^//?" we have that 

2n+2 

^ = l-n^— — . (38) 

SO for R fi^ a. shift in the scalar curvature of order /x^ implies a shift in i/j of order one 
which corresponds to an energy scale of order Mp for the canonically normalized field 
0. In these models, once the background scalar curvature reaches its vacuum value it 
stays locked into this i? ~ /i^ regime, and we would need extremely energetic processes 
involving strong fields (0 > Mp) to get it out of this "locked phase" locally. 

3.2 Fine-tuned case 

In p!T1[T2] several models are proposed for which the second derivative of F is zero on a 
particular low curvature (de Sitter) background. It is argued that since the mass of the 
scalar ~ 1/F"{R) now goes to infinity, its action range goes to zero and the model 
becomes compatible with the Solar System experiments. In fact, to be compatible with 
the current limits from fifth force search experiments, rus > 10~^eV would already be 
fine [H]. For instance in the case (|35|) one has 

2 (i?o-3a/i2) ^ ' 

and one finds that > lO^'^eV for 

i?o = 3a/i^(l ± 10-*'°) . (40) 

This already shows the problem with this proposal. If the background curvature devi- 
ates only minutely from this fine-tuned value, the scalar would become light again, in 
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conflict with experiment, as we already noticed in |15]. However, one might still think 
that for these particular fine tuned backgrounds the static solution corresponding to a 
spherically symmetric mass distribution would be very close to the GR solution, cor- 
responding to (129|) with large m^. But, just as in the previous subsection, we should 
check the weak field conditions to make sure that we can actually use the linearized 
solution (129|) . Let us therefore consider the form of Vp{ip) in the neighbourhood of the 
background for which F" = (and where —>■ oo). We will denote by Rt and ipt the 
values of R and ip in such background. As we said, even on this fine-tuned background 
the condition fl33|) will still be satisfied for some n > 2. For instance in the cases fl35|36l) 
this condition is satisfied for all n > 2. In particular we have F'"{Rt) ~ l//i^ and one 
can show from fl6H71) and ffTTj) that: 

^ F"'{R,)^-^-^, (41) 
Vp(V^) - Vp(^t) ^ -/'(7^(^i))(^W-^(^^)) 

- ±v/2^.^l^v^ (42) 



\F"'iRt)/ 

where the approximations are valid for (R — Rt) <^ /i^ or equivalently (ip — tpt) <^ 1 and 
the sign in (H2!) depends on the chosen branch for TZ. Looking at the approximate form 
of Vp for ip close to tpt it is obvious why the Taylor expansion over a background with 
ipo = V't breaks down, since Vp is non-analytic at this point. This should not come as a 
surprise, because we know that when F" becomes zero one can not invert the relation 
f'{R) = ip and TZ gets a branch point. In any case, one can see that the weak field 
conditions fl30|) for the expansion on a background with curvature Rq close to Rt now 
simply reduce to the condition for the validity of the expansion of Vp{ipo{l + in 
powers of ip, which are 

^-^^o-A-T^jj^TTT^TTT^-l^X -(—X . (43) 



F"'{Rt)V^{^o)^ \mj \ms 

Here rUg corresponds to the mass of the scalar in the background with ip = ipQ. Plug- 
ging here the expressions for the linearized solutions for ip we find that this condition 
becomes 

^^<^. (44) 
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for objects with size r0 -C ^ , and 

G;vP < ^ , (45) 

for objects with size ^ ^7^- And it is apparent that if we take > 10~^eV 
we can not use the hnearized solution in any situation. This failure of the weak field 
expansion is also manifested in the very low strong coupling scale that we now find: 

A,~M,(|)\ (46) 

which is even smaller than the Hubble scale (< H^/Mp) for > 10~^eV. 

To find the true static solution corresponding to a spherically symmetric mass 
distribution with density p we should look at the full equation ffTOj) . Let's take for 
instance the specific background with curvature Rq = Rt- If this static background is 
solution of the EOM we have that Vp{ipt) = 0, so using the approximate expression for 
Vp, eg. (1421) . the scalar equation (ITOl) now becomes: 

±(^^Y = -^-.G^^f~p. (47) 

Notice that this is valid as long as ip, h^^, <^ 1. The correct approximate form of the 
equation does not correspond to a Klein-Gordon equation with a large mass term, 
but an unusual power ip^^"^ of the fluctuation appears in the potential term. Had we 
inappropriately used a Taylor expansion of the potential we would have erroneously 
concluded that the mass term appearing in the equation diverges. Imagine now that we 
have a solution of this equation that matches the asymptotic background at infinity, 
so — s> for r — s> oo and corresponds to the spherical source with mass density p. 
Keeping the correct approximate form of the equation we can see that the Laplacian 
will dominate over the potential term at short distances as usual and the solution close 
enough to the source will be the one of a massless scalar to a good approximation: 

^-^^C. ,48) 

We have not fixed the constant C because its precise value would depend on the 
non-linear interactions of ip that will be important at long distances. By comparing 
the terms in the Laplacian that we have solved for, with the potential term that we 
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have ignored, we see that the solution corresponding to a massless scalar will be a 
good approximation as long as r ^ (G^M/ n^Y^^ if C < (hGnM)^^^ or as long as 
r < (G'ArM/(VC'/i^))^/^ if C > {iiGnMY^^. So we see that regardless of the precise 
value of C (that, as we said, should be determined by matching with the asymptotic 
background), the solution to the full equation fjTr|) in this fine-tuned case is the same 
as the one we found in the generic case for small distances. 

Another specific problem with models of the type fl36|) . that include terms like ~ B? 
in the action, is that one does not even recover Einstein gravity for high curvature 
backgrounds. One can see this immediately from the Lagrangian because it can be 
approximated by / ~ R+aR^ / /i^ when R ^ /i^. In terms of the extra scalar fluctuation 
this means ml ~ Rq for large background curvatures Rq ^ /i^. Notice that now 
A, ~ M^ff, so that we can indeed trust the linearization, and we would find an extra 
long range force when performing local experiments. In the next section we will show 
how the extra scalar fluctuation disappears for large background curvatures, if the 
action does approximate the Einstein-Hilbert action, f — R for R^ fi^. 

Finally, we have only explicitly considered the case with F"'{Rq) equal to its 'nat- 
ural' value. One could imagine an even more fine-tuned case where F^'^^Rq) < H^-^"^ 
for all n < m and (133!) is satisfied only for n > m. In that case one finds for ip — ipt ^ 1: 



v,(^) - VM ~ -"'(^ - i^tY^^"'-'^ , (49) 

resulting in a strong coupling scale 



'Ho 



2 



2m — 2 

m-2 



A, ~ ^ , (50) 



so we see that in any case A^ < Mp (j^^ (~ Hq for iris ~ 10 ^eV). The correct 
approximate equation for the scalar, eq. (H7|) . would now generalize to 

+'-( ^^-^y-^A - .^G^ff, (51) 

V ^ + 3 F(«)(i?0 ) ~ 3""^^ ^' ^^^^ 

and the same argument can be used to show that the solutions will approach those of 
the scalar-tensor theory with a massless scalar at short distances. 
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4 The spherically symmetric solutions on high cur- 
vature backgrounds 

4.1 The Chameleon effect 

Let us begin this section by briefly reviewing the Chameleon mechanism for hiding 
the effects of a field that is otherwise very light in vacuum. Usually one assigns a range 
to the force mediated by a given field according to its mass, because when r > the 
potential produced by the source gets an exponential Yukawa suppression. However, 
when non-linear interactions are important there are more possibilities. Imagine that 
we have a scalar field with an arbitrary potential in vacuum, V{(f)), and a coupling to 
the trace of the matter EMT likqj AC = a{(j))T = —a{(j))p. For finding static solutions 
we should solve the equation 

VV = a'(0)p + ^'(0). (52) 

Notice that the solution of this equation is unique for a given source and asymptotic 
boundary conditions. To understand the Chameleon effect we will consider two different 
spherically symmetric situations: in one p = everywhere except in a small region 
r < Tq where it is constant, while in the other p is a constant everywhere. 

Let's deal now with the first case. If we can use a weak field expansion the solution 
will be completely analogous to the solution of the linearized EOM of f{R) gravity 
presented in the second section but let us briefly repeat it here. The asymptotic value 
of the field, (po, is such that it minimizes its vacuum potential, V'{(j)Q) = 0. Then in 
the region where p ^ the field finds itself in a non-equilibrium position, and it will 
acquire a non-trivial profile. In a weak field expansion in powers of the fluctuation 
= — 00) the linearized equation becomes 

VV=^ + m^(0-0o), (53) 

where = V"{(f)o) and P = MpQ;'(0o). The solution outside the source (r > r©) is 

p-rrtsr 

0(r)=0o + Ci , (54) 

r 

where we have taken into account the asymptotic boundary conditions and Ci is a 



< for all realistic cases, so it is convenient to define p = —T > 0. 
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constant to be determined. The solution inside the source is 

f3p 



^^'■) = - ^rj^ + + C3 . (55) 

MpTTig r r 

Now we can determine the constants Ci by imposing </)'(0) = and continuity of (f) and 
its first derivative at r = r©. Doing this, in the hmit r© <^ we get the usual result 

c = = r56) 

^ 3Mp AwMp' ^ ' 

where M = ^irr^p. We see how the localized energy-momentum density sources the 
field and makes it acquire a non-trivial profile outside the source. 

However, in the second situation where p is a constant everywhere, this term can 
be seen as another contribution to the potential of (p. So the field will not acquire 
a non-trivial profile but will simply set to the minimum of its p-dependent effective 
potential Ve// = V{(f)) + a{4>)p- So the solution will be = 0^ where (ps satisfies 
V'{<j)s) + a'{<j)s)p = 0. 

The Chameleon effect will take place whenever the second situation is a good ap- 
proximation for the solution inside a localized spherically symmetric source. When this 
is the case, the region inside the source where the field is settled to the minimum of its 
effective potential will not source the field outside. The "thin shell" will be the region 
near the surface of the massive body where the field has a profile interpolating between 
its equilibrium positions, 0o ^^nd 0^. If this region has a thickness given by Ar, the 
effects of the force mediated by this field will be suppressed as long as Ar/r© <^ 1. For 
the thickness of this "thin shell" one finds the approximate expression [5]: 

Ar (00 - 0s) 



r© QP^^Mp' 



(57) 

where = G^M/vq is the Newtonian potential at the surface of the body. So the 
condition for the Chameleon mechanism to hold will be that the difference of the values 
of the field that minimize the effective potential inside and outside the source should be 
much smaller than jS^NMp. When non-linearities are negligible we can approximate 
the potential and the coupling function by = m^0^ and a = j3o(j)/Mp. It is easy 
then to compute 0^ and applying the previous formula we see that there will be a thin 
shell only when r© ^ itl~^, and the effects of the field will be hidden only for distances 
larger than the inverse mass, as expected. But when non-linearities are relevant we can 



17 



have a Chameleon effect even if Tq <^ rnj^. Notice that this necessarily implies the 
breakdown of the weak field expansion. In this case only the mass contained within the 
thin shell will contribute to the field outside the source so we can estimate the solution 
for the scalar field outside the source as [9] 

'^-^^'T"^ -^^ ' 

47rMp r Tq AnMp r 

where Mf^ is the mass contained within the thin shell. Through this mechanism a field 
that is very light in vacuum could have passed experimentally undetected. 

4.2 The Chameleon effect in f{R) gravity 

We are now in a position to assess what are the necessary properties for an f{R) action 
in order to get a Chameleon effect for the extra scalar degree of freedom. To parallel the 
Chameleon literature [9] we will work with the canonically normalized field defined 
in eq. (1). For f{R) theories one can see from the action (1121) that the coupling of the 



scalar field to matter is given by AC = e"^'"^T/4 = -e~^'"^p/4 where /t"^ = ■y/3/2Mp 
and T is again the trace of the EMT. Qualitatively, what we need to get a Chameleon 
effect is that as p becomes large, the value of that minimizes the effective potential, 
Veff = V{(j)) + e~'^'^'^p/4:, depends very weakly on p. If this is the case the difference 
00 ~ 0s in the estimation of the thin shell thickness of the previous section will be 
small when the asymptotic background energy density is large and massive bodies will 
indeed develop a thin shell, so that the extra force becomes negligible. When this is 
the case the scalar curvature will follow roughly the Einstein equations, R ~ p/M^. 
The scalar curvature is however given in terms of trough its definition eq.([T]), so 
the requirement that as p becomes large the variation of the equilibrium value for 
becomes small can be rephrased as the requirement that as R becomes large, f'{R) 
becomes roughly constant. So, for modifications of the GR action characterized by a 
curvature scale /i^, we will have a Chameleon effect at high values of the background 
energy density, po ^ A^^^p, (or curvature Rq^ p?) if 

f\R) ~ 1 when R > /il (59) 

This qualitative discussion agrees with the naive expectations that one could have, 
since we are simply saying that the effects of the extra degree of freedom will be 
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hidden when the background curvature is large if the action looks like the Einstein- 
Hilbert action when the curvature is large. However, in order to consider these large 
curvature backgrounds we need them to be relatively stable. This will give a condition 
on the departure from the Einstein-Hilbert action for large curvatures. From (159!) we 
see that there are essentially two possibilities. Either f"{R) will be close to zero and 
positive; or f"{R) will be close to zero but negative. In the latter case we find a large 
tachyonic instability for the scalar fluctuation (m^ ~ l/f"{R)), resulting in a decay 
of the large curvature background. For instance, the models f{R) = 1 — ^^^^^ / R^ 
that were proposed originally [16], have ~ —R{R/fi^)"''^^ for large curvatures, so 
for such models those backgrounds are unstable. For this reason they fail to give a 
realistic expansion of the Universe at early times, as noticed in |5]. This also agrees 
with the results of [T7], where it was shown that FRW solutions for these models 
never attain large curvatures, clearly in conflict with a conventional matter dominated 
expansion. 

The story is completely different for the models with positive values of f"{R) at 
large curvatures, that we will consider from now on. The scalar fluctuation now has a 
large positive mass squared and the backgrounds are stable. So we can safely assume 
that the FRW solutions will only start to differ from those of Einstein gravity at 
the current epoch, when if ~ i^o ~ One might think that the recovery of Einstein 
gravity for large background curvatures now simply happens because the mass becomes 
large, reducing the range of the extra force. However, this assumes that the weak field 
expansion remains valid. But, as we preempted in the beginning of this subsection, the 
weak field expansion will in fact break down for realistic situations, and the recovery 
of Einstein gravity happens through the non-linear Chameleon mechanism. For a 
massive body of mass M, with radius Tq and density p = 3M/(47rrQ) we then get the 
following picture. For low values of the asymptotic background energy density po the 
gravitational field of the body will be of the scalar-tensor type, with an extra force 
with long range l/m^. If we now increase the background density and curvature, the 
range of the extra force will decrease. But then at some point, when l/m^ is still 
larger than Tq, the linearized solution breaks down and the body develops a thin shell, 
suppressing the extra force. Let us now illustrate this for some particular form of the 
function /. We will consider functions such that when R ^ fi^ can be approximated 
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b}|^ f{R) — R + /i^"^^/-R". The relation of the scalar field with the curvature will then 
be eiven by 

,,2n+2 

e"*-l-^Tr^- (60) 

Since we are assuming that R is positive and much bigger than fi we will be interested 
in the behavior of the potential for negative values of close to zero. Our assumptions 
imply then that in this region the effective potential for the field can be approximated 
by 



For Po ^ ^pP^ this potential is minimized for 



4e-''^'-^^pV(-^)^- (62) 



K(j)Q ~ — n — - — . (63) 
V Po / 

We see how the equilibrium value of the field simply gets closer to zero as we increase 
the background energy density. As we just said, this will give a large mass for the 
scalar fiuctuations on these large curvature backgrounds: 

(\ n+1 / f? \ 

However, at the same time this gives a low strong coupling scale for the weak field 
expansion. Indeed, the expansion of the potential fl62l) in powers of the fiuctuation cj) 
breaks down for |0| > |0o|, so we get a strong coupling scale 

A, ~ |0o| -mJ^^J , (65) 

and the linearized solution for 0, that we get from eg. (1531) with fi = — /tMp/2, is only 
valid when 

2GnM , , , /M,V^"^' 



«:0(ro) ^ — ^ « ~ . (66) 

Sr© V Po 



^Notice that we are not making any assumption here about the properties of the vacuum solutions 
(where R ^ ji^) since we are just assuming a functional form for / that holds in the high curvature 
limit i? > 
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For a given source, this will be the case for background densities po smaller than a 
certain critical density pc given bjc^ 

where $Ar stands again for the Newtonian potential at the surface of the body. If the 
background density is much larger than pc the linearization breaks down and the full 
non-linear equation is approximately solved by the Chameleon thin shell solution (!58|) . 
Indeed, from eq. fl571) we find that the body develops a thin shell if: 

= - « 1- 68 

\ Po J \PoJ 

Notice that massive bodies with strong gravitational fields in their surface {i.e. large 
$Ar) will develop a thin shell more easily than smaller ones for a given asymptotic 
background curvature. We also see how, when we take the limit p — > leaving po fixed, 
all sources will develop an (infinitely) thin shell and we recover the solutions of GR, in 
which the field (curvature) is a constant outside the source and jumps discontinuously 
(Ar ^ 0) to a different value inside the source. So we see that the true limit to GR in 
these models would be ^ — * 0, where po is a background energy density. 

In these estimations, if we are dealing with the cosmological background, we can 
approximate ~ (1 + 2;)^. Assuming then an adiabatic evolution {i.e the solution is 
always taken to be the equilibrium, static one for the prescribed asymptotic value of the 
curvature,) we can estimate using eq. (1681) at what cosmological time a given isolated 
source would change its gravitational field from GR to scalar-tensor. For instance a 
star has typically $jv ~ 10"^, while a galaxy or galaxy cluster can have $iv in the range 
~ 10^^ — 10^''. It is clear thus from (1551) that we can expect that the gravitational field 
of these astrophysical sources would have changed to the scalar-tensor type already at 
high redshifts, when [1 + z) ~ 

^Q4/(3n+3) _ ^Q7/(3n+3)_ Notlcc that if we had used the 
invalid linearized result fl54ll56p . we would simply assign a range to the force mediated 
by the field given by its inverse mass, ~ and we could estimate the redshift 

dependence of this distance as rs{z) ~ m~^{z) ~ ^^""^(1 + z)^^("+^)/^. So we would 
conclude that the field would have very long range even at high redshifts. However, 



""^"We are assuming here that the density of the body, is p 3> Pc', one can show from ([M)) that this 
indeed imphes that l/ms(pc) S> tq. 
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here we have shown that non-hnear interactions provide a further suppression of the 
effects of the scalar field through the Chameleon effect. 

We should point out here that in this section we have just presented ordcr-of- 
magnitude estimations for the necessary conditions to recover Einstein gravity. But 
since these solutions are non-perturbative, in order to study its behavior at a more 
quantitative level when for instance Ar/r© is not very small, a numerical integration 
of the equations would be mandatory. Also, we have assumed an adiabatic evolution 
in the estimation of the cosmological time at which the gravitational field of a given 
source would change from GR to scalar-tensor. Under this assumption the models 
would be ruled out, since we would be in the scalar-tensor regime at present. But 
to study the time-dependence of these solutions, thereby checking explicitly if their 
evolution is really adiabatic, one might have to resort again to numerical methods. On 
the basis of our analysis we can not exclude that models exist for which the GR to 
scalar-tensor transition is slow and non-adiabatic. In this case such f{R) models could 
perhaps be brought into accord with observations. 

5 Conclusions 

In this paper we have studied the solutions corresponding to spherically symmetric 
sources in f{R) theories of gravity. We have started by clarifying the range of vahdity 
of the linearized expansion, giving the conditions that have to be satisfied for this 
expansion to be valid for an arbitrary function /. Then we have shown that for the 
models that represent a modification of gravitational dynamics only at low curvatures, 
the linearized expansion on vacuum breaks down as long as the scalar curvature departs 
significantly from its vacuum value. These models are characterized by a function that 
can be approximated by the Einstcin-Hilbcrt action (/ ~ R) when R ^ fi^, but is 
non-trivial when i? ~ /x^. However, when we fix our asymptotic boundary conditions 
for the curvature to the vacuum value, in most cases the linearized solutions (that are 
in clear conflict with Solar system experiments) are vahd as long as GnM/t <^ 1, 
the same condition that one finds in GR. Indeed, the strong coupling scale of this 
linearization is the Planck mass, as in GR. This implies that, in these models, once the 
scalar curvature diminishes to a value ~ //^, it gets "locked" into this value everywhere, 
and one would require very strong gravitational fields or Planck scale energies to be 
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able to get out of this phase locally. It is worth to mention here that the situation 
is completely different for models that include inverse powers or logarithms of other 
curvature invariants beyond the scalar curvature. In particular if we include invariants 
that do not vanish in the Schwarzschild solution, even in the case when the linearized 
solutions are the same, the strong coupling scale for the linearization on vacuum is 
significantly smaller, ~ (yU^Mp)^/^, as we showed in [18]. For those models one 
can never use the linearized solutions, and one does recover the GR solutions at short 
distances as required by Solar System or laboratory experiments [18] . 

Also, we have seen why raising the mass of the scalar field in vacuum by adding, 
e.g. an term to the action does not imply compatibility with local tests of gravity, 
as claimed in [12]. In fact, it is easy to see that those models are ruled out because 
when adding positive powers of R to the action we are modifying gravity also at 
high curvatures. For instance in the R^ case one can approximate the action by / = 
R + R^/rn^ with m ~ /i in the high curvature regime {R ^ ^'^). In those situations 
the scalar field becomes effectively massless and its effects would have been observed. 
Furthermore, we have argued that even the solutions in vacuum are of scalar-tensor 
type for observationally relevant distances in these cases. 

On the other hand, for the models that modify gravity only at small curvatures, 
when the asymptotic value of the curvature is large {i.e. ^ /i^ because there is a large, 
constant background energy-momentum density), the linearization breaks down and 
the behavior of the extra scalar field, that we can associate to the extra degree of 
freedom contained in the metric, is governed by non-linear dynamics. This observation 
allows one to recover quantitatively what intuitively could seem obvious: the solutions 
of the theory approach those of GR when i? ^ /i^. This is achieved through the 
so-called Chameleon mechanism [Hj: whenever there is a localized massive body in a 
high curvature background, the scalar field quickly adopts a new equilibrium position 
inside it as a consequence of its non-linear interactions even if its mass is small on the 
background. As a consequence only a "thin shell" of matter in the surface of the body 
acts as a source for the field outside. This "thin shell" is the region near the surface 
where the scalar field interpolates between its equilibrium positions. We have estimated 
this thickness for some forms of the function /, giving a quantitative estimation of the 
necessary value of the background curvature in order for this effect to take place for 
a body of given mass and radius. Using this estimation, and under the assumption 
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of an adiabatic evolution, we have argued that in the apphcation of these theories to 
the Universe one should observe a transition from GR to scalar-tensor gravity as the 
Universe expands, taking place already at high redshifts for most astrophysical sources. 
However we have not ruled out the possibility that this gravitational decay could be 
slow and non-adiabatic for some models. If the Earth and its environment would 
still remain in the GR regime in some cases, those models would be compatible with 
local tests of gravity. But a quantitative determination of the time scale associated 
to the decay of the scalar curvature to its equilibrium value ~ /i^ inside sources as 
the asymptotic background curvature approaches its vacuum value would require the 
study of time-dependent non-perturbative solutions which lie beyond the scope of the 
present paper. 

Finally, our analysis also allows us to make some more general statements on the 
applicability of the linearization for these theories. For low curvatures of the back- 
ground -Ro ~ one can safely use the weak field expansion in the same situations 
as for GR. However, for large curvatures, one has to be careful when using this ex- 
pansion. For instance, the linearized equations for the cosmological perturbations that 
were used in [5] will become invalid at early times, at least if one considers realistic 
{i.e high curvature) FRW backgrounds. The recovery of Einstein gravity is essentially 
non-perturbative as we illustrated in the previous section. And in such models we can 
expect that in general there will be a range of background curvatures for which the 
linear solutions are those of scalar-tensor gravity, whereas the true non-linear solutions 
are in agreement with GR. 

Acknowledgements 

We acknowledge discussions with Eanna Flanagan and Luca Amendola. K.V.A. is sup- 
ported by a postdoctoral grant of the Fund for Scientific Research-Flanders (Belgium). 

References 

[1] B. Whitt, Phys. Lett. B 145 (1984) 176; J. D. Barrow and S. Cotsakis, Phys. Lett. 
B 214 (1988) 515; K. i. Maeda, Phys. Rev. D 39 (1989) 3159; D. Wands, Class. 
Quant. Grav. 11 (1994) 269 |arXiv:gr-q c/9307034]. 



24 



[2] T. Chiba, Phys. Lett. B 575 (2003) 1 |arXiv:astro-ph/0307338j . 

[3] D. N. Vollick, Phys. Rev. D 68 (2003) 063510 |arXiv:astro-ph/0306630|; 



E. E. Flanagan, Phys. Rev. Lett. 92 (2004) 071101 |arXiv:astro-ph/030811l| . 
[4] S. Capozziello, V. F. Cardone and A. Troisi, Phys. Rev. D 71 (2005) 043503 



arXiv:astro-ph/0501426|; S. Nojiri and S. D. Odintsov, Phys. Rev. D 74 (2006) 



086005 |arXiv:hep-th/0608008| . 
[5] Y. S. Song, W. Hu and I. Sawicki, |arXiv:astro-^h706 10532] 



[6] G. J. Olmo, Phys. Rev. Lett. 95 (2005) 261102 |arXiv:gr-qc/0505101| , Phys. Rev. 
D 72 (2005) 083505. 



[7] A. L. Erickcek, T. L. Smith and M. Kamionkowski, arXiv:astro-ph/0610483 



[8] J. A. R. Cembranos, Phys. Rev. D 73 (2006) 064029 |arXiv:g r-qc/0507039|; 
C. G. Shao, R. G. Cai, B. Wang and R. K. Su, Phys. Lett. B 633 (2006) 164 



[arXiv:gr-qc/0511034 



[9] J. Khoury and A. Weltman, Phys. Rev. Lett. 93 (2004) 



171104 |arXiv:astro-ph /0309300|, Phys. Rev. D 69 (2004) 044026 
[arXi v:astro-ph/03 09411|; P. Brax, C. van de Bruck, A. C. Davis, J. Khoury and 
A. Weltman, Phys. Rev. D 70 (2004) 123518 far Xiv:astro-ph/0408415| ; D. F. Mota 
and D. J. Shaw, Phys. Rev. Lett. 97 (2006) 151102 [arXiv:hep-ph/0606204| , 



arXiv:hep-ph/0608078 



[10] A. D. Dolgov and M. Kawasaki, Phys. Lett. B 573 (2003) 1 
[arXiv:astro-ph/0307285 |; V. Faraoni, [a?Xiv:astro-ph/0610734 



[11] R. Dick, Gen. Rel. Grav. 36 (2004) 217 [arXiv:gr-qc/0307052l . 

[12] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68 (2003) 123512 
[arXiv:hep-th7030728 8|, Gen. Rel. Grav. 36 (2004) 1765 | arXiv:hep-th/0308176 . 



[13] B. Bertotti, L. less and P. Tortora, Nature 425 (2003) 374. 

[14] E. G. Adelberger, B. R. Heckel and A. E. Nelson, Ann. Rev. Nucl. Part. Sci. 



53 (2003) 77 |arXiv:hep-ph/0307284j ; E. G. Adelberger, B. R. Heckel, S. Hoedl, 



C. D. Hoyle and D. J. Kapner, arXiv:hep-ph/0611223 



25 



[15] I. Navarro and K. Van Acoleyen, JCAP 0609 (2006) 006 |arXiv:gr-qc/0512109| . 

[16] S. M. Carroll, V. Duvvuri, M. Trodden and M. S. Turner, Phys. Rev. D 70 (2004) 
043528 |arXiv:astro-pli/0306438 1 . 



[17] L. Amendola, D. Polarski and S. Tsujikawa, arXiv:astro-ph/0603703 



[18] I. Navarro and K. Van Acoleyen, Phys. Lett. B 622 (2005) 1 |arXiv:gr-qc70506096] , 
JCAP 0603 (2006) 008 [arXiv:gr-qc/0511045 



26 



